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Abstract. We present a detailed model of an extreme slim disk (ESD) with a hot corona around a massive black hole with 
dimensionless accretion rate rh in the range 2.5 <C rh < 100. We assume that a fraction / of the gravitational energy is 
released in the hot corona and the rest is released in the ESD. The energy equation of the ESD is dominated by advection and 
the spectrum shows a broad "hump" caused by saturated Comptonization with monochromatic luminosity given by vL v oc u°. 
This relationship enables us to estimate the black hole mass from the ESD luminosity. The spectrum of the hot corona is 
sensitive to the parameter fM and the viscosity and shows a Comptonized power-law with a high-energy cutoff. The model is 
used to explain the spectral properties of narrow line Seyfert 1 galaxies (NLSls). In particular, it can explain the spectrum of 
extreme NLSls like NLS1 RE J 1034+396. Our spectral estimate of the black hole mass in this source is in good agreement with 
the mass obtained by applying the reververation mapping correlation. We also examine the Eddington ratios in a large NLS 1 
sample and find that most objects show super-Eddington accretion rates. We argue that soft X-ray humps in NLSls are natural 
consequences of super-Eddington accretion in such objects. 
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1. Introduction 

Narrow line Seyfert 1 galaxies (NLS1) receive much atten- 
tion because of their unusual properties: the H/3 line is nar- 
row (FWHM H/3 < 2000 km/s) and relatively weak (less 
than a third of the intensity of the [OIII]A5007 line), and the 
optical Fell lines are very strong. These characteristics lo- 
cate NLSls near the extreme value of the Boroson and Green 
(1992) Eigenvector 1. Soft X-ray excess and the tendency of 
the spectrum to flatten at low energy, have been found in sev- 
eral NLSls, e.g. PKS 0558-504 (O'Brien et al. 2001), Mrk 
766 (Boiler et al. 2001), PG 1244+026 (Fiore et al. 1998), PG 
1404+226, PG 1440+356, PG 1211+143 (George et al. 2000), 
IRAS 13224-3809 (Vaughan et al. 1999), 1H 0707-495 (Boiler 
et al. 2002, Dewangan et al. 2002), RX J1702.5+3247 (Gliozzi 
et al. 2001) and Ton S180 (Turner et al. 2001a, Turner et al. 
2002). BeppoSAX and ASCA observations of Akn 564 and RE 
J 1034+396 clearly show these extreme properties including a 
soft X-ray hump with a very flat spectrum (vF v = const., 
Comastri et al. 2001; Puchnarewicz et al. 2001; Turner et al. 
2001b). Other unusual X-ray properties are discussed in Brandt 
et al. (1998), Comastri et al. (2001), Pounds, Done & Osborne 



Send offprint requests to: J.-M. Wang, 

e-mail: wang@astro . uni-tuebingen . de 



(1995) and Collinge et al. (2001). Earlier suggestions attribut- 
ing the soft excess to blends of emission lines are not supported 
by recent ASCA , XMM and Chandra observations (e.g. Turner, 
George & Netzer 1999; Puchnarewicz et al. 2001, Turner et al. 
2001a, b, Collinge et al. 2001). 

According to some theories, NLSls contain relatively small 
(considering their luminosity) black holes (BHs) with very high 
accretion rates (Laor et al. 1994, 1997, Boiler et al. 1996) ex- 
ceeding, in some cases, the critical Eddington rate (Mineshige 
et al. 2000, Collin et al. 2002). This hypothesis is supported by 
comparing NLS 1 spectra with the spectrum of Galactic black 
hole candidates (Pounds et al 1995), and by very few direct 
measurements of the BH mass in a handful of NLSls (e.g. 
Peterson et al., 2000). The observed properties of NLSls, in 
particular the soft X-ray excess and the small BH mass, suggest 
that they may produce much of their emission due to processes 
in high accretion-rate disks with hot coronae. It is therefore im- 
portant to calculate models of such systems. 

Thin accretion disks are thought to produce the optical-UV 
continuum radiation in Broad Line Seyfert 1 galaxies (BLSls; 
see Shields 1978; Malkan & Sargent; 1982, Malkan 1983; Laor 
and Netzer, 1989). Indeed, thin disk spectra, based on the stan- 
dard disk model (Shakura & Sunyeav 1973) with Eddington ra- 
tio L/L^dd < 0.3, give reasonable fits to the spectrum of many 
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Table 1. Model List of Emergent Spectrum from Accretion Disk with High Accretion Rate 



reference 


disk struc. 


radiation 


Comp. 


z— structure 


GR 


corona 


note 


CT95 


global 


BB, f-f 


Yes 


No. 


No. 


No. 


shocks (?) 


SAM96 


global 


MBB 


No. 


No. 


No. 


No. 




WSLZ99 


global 


f-f, BB 


No. 


Yes 


No. 


No. 




MKTHOO 


global 


BB 


No. 


No. 


No. 


No. 




WN 


self-similar 


f-f, BB 


Yes 


No. 


No. 


Yes 





CT95: Chakrabarti & Titarchuk (1995); SAM96: Szuszkiewicz, Abramowicz & Malkan (1996); 
WSLZ99: Wang et al. (1999); MKTHOO: Mineshige et al. (2000); WN - the present work 
f-f: free-free emission, BB: black body radiation, MBB: modified black body radiation by elec- 
tron scattering, Comp.: Comptonization; GR:general relativity. 



BLSls (Czerny & Elvis 1987, Wandel & Petrosian 1988, Sun & 
Malkan 1989; Laor & Netzer 1989; Laor 1990; Ross, Fabian & 
Mineshige 1993, Storzer, Hauschildt & Allard 1994; Shimura 
& Takahara 1994, Dorrer et al. 1996; Sincell & Krolik 1997; 
Hubeny & Hubeny 1997; Kawaguchi et al. 2001). However, 
there were very few attempts to calculate the spectrum of high 
accretion-rate disks which we call either "slim disks" (SDs, af- 
ter Abramowicz et al. 1988) or "extreme slim disks" (ESDs), 
depending on their accretion rate. The few models that are 
available in the literature are shown in Table 1 where we list 
their more important ingredients. Some common properties re- 
garding the emergent spectrum are the flattening at soft X-ray 
energies, the high energy cutoff which is almost independent 
of the accretion rate, and the fact that the total luminosity is 
saturated (Wang & Zhou 1999; Wang et al. 1999; Fukue 2000; 
Mineshige et al. 2000). 

Another ingredient in some disk models is a hot corona. If 
exists, such a corona must be complex, as is evident from the 
highly variable X-ray spectrum. The geometry and formation 
of the hypothetical corona are poorly understood (see a review 
by Collin 2001). There are several published models of stan- 
dard, thin accretion disks with hot coronae (see Shapiro et al., 
1976; Liang & Price, 1977, Svensson & Zdzarski 1994, Janiuk 
& Czerny 2000), but, so far, no theoretical model of a high 
accretion-rate disk with a hot corona. 

In this paper we study ESDs with hot coronae in order to 
explain the X-ray spectrum of NLS Is. The paper is arranged as 
follows: §2 gives the basic equations and explains the method 
of calculations. §3 discusses the spectral properties of such sys- 
tems and in §4 we examine the Eddington ratio in several ob- 
served systems and apply our model to the soft hump observed 
in RE J1034+396. Finally, in §5 we discuss implications to fu- 
ture models and observations of NLS Is. 



2. The structure of ESDs with hot coronae 

2. 1. Extreme slim disks 

The standard accretion disk model (Shakura & Sunyeav 1973) 
breaks down when the accretion rate approaches the Eddington 
rate. At this limit, the disk becomes "slim" and radial advection 
cannot be neglected any more. A large range of accretion rates 
(M) and relative rates (m) of such systems has been considered 
(Abramowicz et al. 1988; Szuszkiewicz et al. 1996, Wang et 



al. 1999, Mineshige et al. 2000; Watarai & Mineshige 2001; 
Ohsuga et al. 2002). Here 



m = 



M 

Mcr 



647rGM B H"i P 



(JtC 



(1) 



or is the Thomson cross section, m p the mass of the proton, c 
the speed of light, Mbh the BH mass and we have assumed an 
accretion conversion efficiency of r/ = 1/16. ESDs are taken 
here to be those cases with 2.5 <C m < 100 (see discussions 
below). In such disks, the emitted photons are trapped in the 
gas due to the large Thomson depth which causes the diffu- 
sion time to exceed the radial motion time (Begelman & Meier 
1982; Wang & Zhou 1999; Spruit 2000; see Ohsuga et al. 2002 
for a more sophisticated treatments of photon trapping). In such 
cases, the surface cooling is less efficient with important con- 
sequences to the emergent spectrum. Below we consider ESDs 
where advection cooling dominates over other surface cooling 
due to such trapping. 

In this work we consider also a hot corona that modifies 
the emergent disk spectrum by Compton scattering, and con- 
tributes to the high energy radiation. Previously published disk- 
corona models were based on the following assumptions: 1) a 
fraction / of the gravitational energy is released in a hot corona 
and (1 — /) is deposited in the cold disk (Haardt & Maraschi 
1991). 2) the vertical equilibrium holds in both the cold disk 
and the hot corona and is described by the vertically averaged 
equations (Svensson & Zdzarski 1994, Janiuk & Czerny 2000). 
3) the hot corona is cooled by Comptonization of soft photons 
emitted by the cold disk including the reflected photons. With 
these assumptions, the structure and the spectrum of slim disks 
with hot coronae can be obtained with the parameter /. There is 
no way to accurately find the value of / which depends on the 
connection between the disk and the corona. We follow a sim- 
ple prescription and treat / as a free parameter of the model. 
The accretion rate in the ESD is thus (1 — f)M. There are 
other constraints on our model, such as the requirement of an 
optically thin corona (see below). 

The basic ESD equations are similar to the SD equations 
for a steady disk in a Newtonian potential of azimuthal symme- 
try and Keplerian angular velocity Ok- We take the kinematic 
coefficient of the shear viscosity to be v = aa s Hd, where a s 
is the sound speed, a the viscosity parameter, and the half 
thickness of the disk. The basic equations are those proposed 
by Muchotrzeb & Paczynski (1982), Matsumoto et al. (1984) 
and Abramowicz et al. (1988), except that we neglect the coef- 
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ficients Bi (which are of order unity) and the accretion rate is 
replaced by (1 — f)M. This gives 



(1 - f)M = 4irRH d p d vf, 



(2) 



Pa 
Pa 



- #d^k> 



(1 - f)M(l - / in ) = ATiR 2 H d aP A , 

^-(tf-nl)R + v? d ^ + ^ d -^ = o, 

p d dR K h> dR p d dR 



(5) 



where Pd is the total pressure, p d the mass density, vf the ra- 
dial velocity of the flow, fl the angular velocity, I the specific 
angular momentum, and l[ n the eigen value of angular momen- 
tum at the inner boundary. The last term in equation (5) is the 
correction for the decrease of the radial component of the grav- 
itational force away from the equator (Matsumoto et al. 1984). 

While the above equations are suitable to describe the slim 
disk-corona coupling, such calculations are beyond the scope 
of the present paper. We only consider Thomson thin coronae 
in the limit /« 1. Despite the small /, the hot corona can be 
a strong emitter since JM is not necessarily small. 

The energy conservation equations should be added to the 
above equations. Advection cooling, which can dominate over 
the diffusion cooling, is given by 

(l-f)M dS (l-f)M fP A \ 

0adv - - Td —^R- 1r = [y d ) L (6) 

(Abramowicz et al. 1988) where S is the entropy, 

£=(4-3/3)7„-(12-10.5/?)7 T , (7) 

7 p = d\np d /d\nR,-f T — d\nT d /d\nR and (3 is the ratio of 
gas to total pressure. In this paper we take (3 — since radiation 
pressure dominates over gas pressure. Thus, £ = 47 p — 12j T . 
For intermediate values of /, the energy equation is Q a d v + 
Qcooi = Qvis' where Q^ is is the energy generation via viscosity 
in the SD and Q^ ool is the cooling from the surface of the SD. 
For ESDs, the condition / « 1 guarantees that advection is 
the dominant cooling process (Q a dv > Qcooi) an ^ practically 



balances Q^ is , i.e., 

Oadv = QL = (1 - /)Qv 



(8) 



where Q v i s is the energy generated by viscosity per unit area 
M(Z-Z in ) / dQ\ nJMfl 2 



Q 



4irR 



dR 



47T 



(9) 



where n = — dlnCl/dlnR and J = 1 — Strictly 
speaking, the radial motion and energy conservation equa- 
tions can be solved given the appropriate boundary conditions 
(Abramowicz et al. 1988). Denoting j p = d In P d jd In R, and 
"f v = d In vf/d In R, the radial motion equation (5) is reduced 
to an algebraic equation 



>-0^-(« 2 -^ 2 +7,K d ) 2 =O. 



(10) 



Neglecting the boundary condition, i.e., setting J — 1 (or 
/in ~ 0), and using eqns. 2, 3, 4, 8 and 10, we get, after some 
algebraic manipulations, the following set of self-similar solu- 
tions: 



(3) P d = 



(1 - QMQk 
AiraR 



Pa = 



7o 2 (1~/)M 
47ra n k R 3 ' 



(4) vf = a 7o 1 RQ k , H d =j 1 R, 



n = 



7o 



a s = 



Rtt k 

7o 



(11) 
(12) 
(13) 



where 70 = (5 + a 2 /2) 1 / 2 is a weak function of the viscosity 
a which is roughly \/5 for a -C 1. It is interesting to note that 
vf, Hd, fi, and a s are independent of /. When / goes to zero, 
the solution reduces to the one given in Wang & Zhou (1999). 

The usual treatment of the inner edge of the accretion disk 
assumes a certain last stable orbit, R ms , depending on the BH 
properties (e.g. R ms = 3R S for a Schwarzschild black hole 
with R s = 2GMbh/c 2 ). This is also taken to be the innermost 
emitting radius i?; n . Most of the calculations presented here as- 
sume a Schwarzschild BH with the corresponding inner radius. 
However, Mineshige et al. (2000) and Watarai & Mineshige 
(2001) have shown that a substantial amount of radiation can be 
produced from inside this value of R ms if the accretion rate is 
larger than 5M C . These authors take R in to be a free parameter 
of the model. Moreover, it has been argued (Krolik & Hawley 
2002) that in real disks, i?; n < R ms , i.e. a substantial fraction 
of the observed radiation is emitted inside the last stable orbit. 
Our calculations cannot include the Krolik and Hawley effects 
and we have adopted the general approach of a free i?; n as in 
Mineshige et al. (2000). We note that the smaller R in is, the 
better is the self-similar approximation since, as explain, for 
this condition to hold we require J ~ 1. 

While surface cooling in ESDs is very small compared with 
the advected energy, the diffuse photons play an important role 
in cooling the corona. We can calculate Q^ ool via the diffu- 
sive approximation. The temperature at the central plane of 
an ESD is given by T d = (SPd/a) 1 / 4 = 6.69 x 10 7 [(1 - 
/)m/m] 1 / 4 o; _1/ ' 4 r _5 / 8 K, where a is the radiation density 
constant and the effective temperature is T c s ~ Td/r V 4 (Laor 
& Netzer 1989). For Thomson scattering, r os = arpdHd, lead- 
ing to 

T cS = 3.64 x 10V 1/4 ™bhV 1/2k - (14) 

where tobh = M-qh/Mq, and r = R/R s . We thus get an ef- 
fective temperature distribution of the form T cS oc r 



-1/2 



as 



in Wang & Zhou (1999) and Fukue (2000). This radial depen- 
dence is caused mostly by photon trapping. The surface energy 
flux,Qi ol ,is 



cool 



cff 



1.0 x 10 26 7^" 1 mBH'" 



erg s 1 cm 2 (15) 



We will show later that the black body approximation used here 
to calculate the emergent flux is justified because of the satu- 
rated Comptonization in ESDs. 

It is important to note that T c g and Q^ ool are independent 
of the accretion rate and the factor / but depend on the mass of 
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the black hole and very weakly on the viscosity a through 70. 
Thus we can estimate the BH mass directly from the spectrum. 

In ESDs, the time scale of photon diffusion in the z- 
direction is T difi ~ iJ d K os E/c whereas the radial transportation 
time scale of the accreting gas is r R <~ ttR 2 Y>/M (Mineshige 
et al. 2000). The condition tr = r d iff can be used to obtain the 
dimensionless photon trapping radius, inside which the photons 
will be trapped in the gas, 



r trap = 7.2 x 10 2 (1 



f) - 
J ' \ 50 



(16) 



The total ESD luminosity integrated over the photon trapping 
region is obtained from the integration of equations (15) 

"(l-/)ml 



L rad w 8.0 x 10 37 m BH In 



50 



ergss 



(17) 



We see that the total radiated luminosity L rad is only weakly 
dependent on the accretion rate to. Comparing with the thin 
disk critical radius k, 14 (Laor & Netzer 1989), inside which 
half the total luminosity is dissipated, we see that for high ac- 
cretion rates, most of the gravitational energy is released within 
the trapping radius and much of it is advected into the black 
hole. Because of that, the calculations of the ESD spectrum ne- 
glect the region beyond r trap that contributes very little to the 
emergent spectrum. 

A necessary condition for our self-similar solution to hold 
is Qadv S> Qcool- This is equivalent to to S> 0.1 8r. If we 
demand that advection dominates surface cooling at r = r^, 
namely at least half the dissipated energy will be advected, we 
get to 3> 2.5 which, using eqn. 16 gives r trap w 36. As ex- 
plained, such a high accretion rate leads to a significant release 
of gravitational energy inside the last stable orbit. Thus, in our 
model we allow R ln < 3i? s . 

2.2. The hot corona 

We assume that the gas in the optically thin corona rotates with 
Keplerian velocity, and that gas pressure dominates over radi- 
ation pressure since the corona is optically thin. Advection is 
mainly due to protons that are heated close to the virial temper- 
ature while the temperature of the electrons is much lower due 
to the inefficient energy exchanges between electrons and pro- 
tons and due to Comptonization cooling. The assumption that 
vertical equilibrium holds means that the height of the corona 
(H c ) is given by H c = c s /Q k , where c s = (P c /p c ) 1/2 oc Tp /2 
is the sound speed. The cooling in the corona is determined 
by Comptonization, and the seed photons are contributed from 
two sources: 1) direct emission from the ESD and 2) ESD 
reprocessing of photons scattered down by the hot corona. 
Considering the fraction r\ of the coronal flux directed towards 
the ESD, and the energy averaged albedo a, Comptonization 
cooling can be written as 

Q c = A[Q d cool + v (l-a)Q c ] , (18) 

where A — exp(y) — 1 is the Compton amplification factor, 
y = t cs S (1 + 5 e ), and S = 4kT e /m e c 2 . Following Zycki, 
Collin-Souffrin & Czerny (1995), we assume 



where a c is the viscosity parameter in the corona, which may 
be different from the viscosity in the ESD. Q adv is the advected 
energy due to the protons in the corona, 

fMc 2 

Qadv = ^ S ( R ), (20) 

where S(R) = d\nP c /d\nR — 2.5d In T p /d In R is a param- 
eter controlling the corona advection (Janiuk & Czerny 2000) 
that should be obtained self-consistently since it is a function of 
radius (see next section). The energy equation in the hot corona 
can be written as 



Qc + Q adv = /Qv 



(21) 



Qc + Qldv = ^ k a c P c H c , 



(19) 



As explained, in such an optically thin medium, the electron 
and proton temperatures can be quite different due to inefficient 
energy exchange. Thus, the hot corona is a two-temperature 
plasma. The net flux from protons to electrons is given by 

Q c - \.hkm-\T p - T e )(l + ~y> ei p c .ff c , (22) 

where v el = 2.44 x lO 21 /^" 1 ' 5 In A with In A w 20 (Shapiro 
et al. 1976). We take 77 = 0.5 and a = 0.2 throughout this 
paper. The free-free emission is j v ss 6.8 x lO -38 ™ 2 ! 1 ^ - 5 
(we neglect the Gaunt factor). The ratio of f-f emission to 
Comptonization obtained from eq.(22) is 5.6i/ig(T e /T p ) <C 1 
if T e -C T p , where vig = f(Hz) /10 19 . Thus free-free emission 
can be neglected in the two-temperature plasma. The structure 
of the hot corona is determined by viscosity heating, and ad- 
vection and Comptonization cooling. 

An important constraint on the value of / is obtained from 
limits on the Thomson depth of the corona. The number den- 
sity near the last stable orbit (R ms ) is n c = f M / Air R^cnip, 
where we take the radial velocity to be the speed of light and 
the height of the corona to be R ms . The Thomson depth is thus 
t c = n c (7T^ms = ffh/Qflo- This means that / < 6j7oto _1 
or else the corona is Thomson thick (Blandford 1990). This 
constraint does not depends on the disk model. To summa- 
rize, the present model is valid when to 3> 2.5 and / < 
0.0375 (to/10) -1 . 

The structure of ESDs with hot coronae is controlled by six 
parameters: a,rh, mbh, Hni a c and /. The first four describe 
the ESD and the other two the hot corona. It is important to 
note that the density and the height of the corona depend on the 
black hole mass, but the corona structure, i.e. T e , T p , and r es , 
is independent of the black hole mass when using the dimen- 
sionless radius r. This implies that the corona structure will be 
similar in stellar and super-massive black hole accreting sys- 
tems. The corona structure is also only weakly dependent on 
the accretion rate. These may be the reasons why the hard X- 
ray spectrum of Galactic black hole candidates at high state is 
quite similar to that of narrow line Seyfert 1 galaxies containing 
much larger black holes. 

We used an iteration method similar to Chen (1995) and 
Janiuk & Czerny (2000) to calculate advection in the hot 
corona. In this procedure, an initial value of advection factor 
<5(i?) at radius R is assumed. We then obtain the structure of 
the corona at R - AR/2 and R + AR/2 and a new 5(R) which 
is used for the next iteration until the value of S converges. In 
practice, the iteration converges very fast. 
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Fig. 1. The structure of hot corona above ESD. The accretion rates, from bottom to top, are rh = 30, 50, 70, 90, 100 M cr and / = 10~ 3 . In 
Fig-la a = a c = 0.1, in Fig lb a c = 0.05 and in Fig lc a c = 0.01. The advective parameter 5(r) is obtained from self-consistent iterations 
(see text). The hot corona becomes optically thick when rh > 100 (for / = 10~ 3 ). Under such conditions, emission from the disk is blocked 
and no hump is observed. Using the dimensionless radius r(= R/R s ), the structure of the corona is independent of the mass of the black hole. 



Figures la-lc show the structure of a hot corona above an 
ESD for various accretion rates M c — fM and viscosity a. 
From Fig- la., we see that the Thomson depth (r cs ), electron 
temperature (T e ) and proton temperature (T p ) increase with 
M c . If we let the factor / increase, the corona will gradually 
become optically thick and blocks the emission from the ESD. 
This gives a constraint on / (/ <C 1). The temperature of the 
electrons in the corona is a weak function of radius. The op- 
tical depth t cs is a power law in radius, and is sensitive to the 
accretion rate fM. 

We find that the proton's temperature is super-virial, T v j r = 
5.45 x 10 12 r _1 K. This temperature is caused by advection, 
in which the energy dissipated by viscosity is stored as en- 
tropy rather than being radiated. Such a process is enhanced 
by the inefficient energy exchange between electrons and pro- 
tons. Advection affects the corona structure, and is found to be 
a source of heating rather than cooling. We also see that 5(r) 
tends to a constant at larger radii. 

We tested the effect of changing the viscosity on the corona 
structure (Fig. lb-lc). We find that the scattering depth t cs in- 
creases slightly with decreasing a c due to the decrease of radial 
velocity with a c , resulting in the increases of scattering depth. 
For a hot corona with lower viscosity and higher accretion rate, 



we find that the location of maximum electron temperature is 
shifting to very large radii (e.g. log r <~ 1.5). These effects have 
no significant observational consequences. 

3. The spectrum of ESDs with hot coronae 

We have calculated spectra of ESDs with hot coronae over a 
large range of BH mass, accretion rate and relative accretion 
parameter /. Like previous works, we find Comptonization to 
be very significant in such systems. The importance of the pro- 
cess can be estimated by considering the parameter 

y = So (^)' (23) 

where the Thomson opacity is k cs = 0.34, the free-free opacity 
is kh = 1.5 x l0 25 pT- 7 / 2 x- 3 (l - e~ x ), and x = hv/kT. To 
evaluate Y, we use the self-similar solution (eqs. 11-13). The 
results are shown in Figure 2 for rh = 30. We find Y » 1 
for high energy photons (x > 0.2) in a large region inside 
the trapping radius, suggesting saturated Comptonization. For 
higher accretion rate, Y is even larger, and Comptonization is 
stronger. Wang et al. (1999) and Mineshige et al. (2000) con- 
sidered the effects of advection but neglected Comptonization 
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1 1.5 2 

\og(R/R s ) 

Fig. 2. The Comptonization parameter Y for m = 30.0 and a = 0.1. 
The solid lines, from bottom to top are: x = hu/kT=0A, 0.2, 0.3, 0.4, 
0.5, 0.6, 0.7, 0.8, 0.9, 1.0. Comptonization is highly saturated above 
the dotted line that indicates Y = 1.0. 

in their calculations. The present paper includes all these pro- 
cesses. Figure 2 shows the Comptonization parameter Y as a 
function of the normalized radius. 

With high accretion rate, ESDs have unique effective tem- 
perature distribution of the form T e g oc r -1 / 2 and relatively 
high temperature. The emergent spectrum differs significantly 
from the spectrum of thin accretion disks with or without coro- 
nae. We use the method described by Czerny & Elvis (1987) 
(see also Wandel & Petrosian 1988) for the spectral calcula- 
tions. The disk structure is self-similar (eqs. 11-13) and the sur- 
face temperature T s is taken to be T e g from equation (14). The 
ESD spectrum is given by integrating over the photon trapping 
region, 



/■-Rtr 

J v = 7T / 

JRin 



B„(T s )f v RdR, 



(24) 



where the factor f v describes the departure from black body 
radiation, including Comptonization, 



U = fs{v) [l-/ th (i/)] + C. 

Here the factors fg and f t h are given by 



fs 



and 



2 1-e 



1 + (Ktot/Kabs) 



1/2 ' 



/t 



th 



exp 



\n(kT s /hv) 



r 2 ln(l + S c + S 2 ) 



(25) 



(26) 



(27) 



where / t h is the fraction of thermalized photons. 
Kg is the total scattering and free-free 



K tot — 

K cs + K ff is tne total scattering ana tree-tree opacity, 
t* = (K tot K ff ) 1 / 2 S is the effective optical depth, and the 
Comptonization factor C is given in Czerny & Elvis (1987). 

Regarding the hot corona, the emergent spectrum can be 
calculated at each radius from the semi-analytical formula in 
Hua & Titarchuk (1995). We assume that the photons are 
Comptonized locally, i.e. we neglect the Comptonization of 



seed photons from other radii. The integration is performed 
over the trapping radius. 

Figure 3 shows spectra of ESDs with hot coronae for dif- 
ferent accretion rates and black hole mass. The most prominent 
feature is a strong soft X-ray energy band. The luminosity of 
the hump is insensitive to the accretion rate because of photon 
trapping and increases only slightly as ~ lnm (eq.17). This 
results from the superposition of saturated Comptonized spec- 
trum at different radii, satisfying the condition Y ^> 1 (Figure 
2). For such a large Y, the factor ,f t h tends to 1 and C tends 
to a constant (eq. 25). With the effective temperature distribu- 



tion T e g oc r 



-1/2 



we get an emergent spectrum that satisfies 



the condition vL v oc v v (eq. 24) and L v is proportional to 
wbh (see below) . This characteristic shape is in rough agree- 
ment with the no Comptonization case (Wang 1999, Wang et 
al. 1999, Mineshige et al. 2000), but the maximum frequency 
is at higher energy. 

The accretion rate determines the trapping radius which 
determines also the shape of the red side of the soft hump. 
Comparing the cases of m = 10 and rh = 40, we find that 
the width of the soft hump is greatly reduced at smaller rh. The 
behavior can be understood by inserting eq. (16) into (14) to 
get the minimum frequency of the hump 



2.82 x 10 16 7o 1/4 m 



1/4 
BH 



-1/2 



Hz, 



(28) 



where we assume / w 0. The maximum frequency of the 
hump, which is determined by saturated Comptonization, is 
independent of rh and is given by ^ max w 2.7kT e g(r- m )/h, 
namely 



1.19 x 10 18 7 1/4 m 



1/4 
BH 



(t) 



-1/2 



Hz. 



(29) 



The width of the soft hump, from v m - m to f max , depends weakly 
on the accretion rate. This is very different from standard thin 
disk models where the maximum frequency and luminosity 
depends on both the BH mass and the accretion rate. The 
monochromatic luminosity in the flat part of the hump is not 
sensitive to the accretion rate but depends linearly on the mass 
of the black hole. This results in a good method to determine 
a black hole mass for super-Eddington accretion flows (see be- 
low). 

The spectrum of the hot corona is determined by two pa- 
rameters: frh and a c . Several examples are shown in Figure 3. 
The general shape is a power-law continuum with a cutoff. The 
spectrum becomes harder with increasing /. This dependence 
can be understood with the help of Figures la-lc that show that 
the optical depth and temperature of the hot corona increase 
with /. The resulting spectrum can easily be understood by 
the simplest version of Comptonization. The temperature deter- 
mines the cutoff energy whereas the photon index p of the spec- 
trum is determined by the scattering opacity and the tempera- 
ture, p = (7 + 9/4) 1 /2 _ 1.5, W here 7 = 4tt 2 /3(t cs + 2/3) 2 S c 
(Sunyaev & Titarchuk 1980). When the accretion rate M c (or 
/) increases, both r os and S c increase, resulting in the decrease 
of 7. Thus, the emergent spectrum becomes harder as the ac- 
cretion rate increases. 
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15 16 17 18 19 20 



log v (Hz) 

Fig. 3. The emergent spectrum from ESDs with hot coronae for dif- 
ferent accretion rates, black hole masses and /. In all cases a = a c = 
0.1. The spectrum is composed of two components: the ESD contri- 
bution at low energy and the hot corona emission at high energy. The 
relative flux of the two depends on the value of / and m which are 
marked on the diagram. Note the linear scaling of the hump luminos- 
ity With 7ttBH. 

The overall spectrum shows a "knee", a transiting from soft 
to hard X-rays. The flux at the knee frequency is contributed 
from the superposition of the Wien tail of the soft hump and 
the Comptonization of the soft hump. The knee's flux increases 
with / and with decreasing viscosity and is correlated with the 
hard x-ray spectral index (Fig. 3). This is caused by the fact 
that the radiation from the hot corona increases with /. Finally, 
The cutoff energy of the hot corona spectrum depends on the 
maximum temperature of the electrons. This cutoff, which is 
below 1 MeV, may be observed by future missions such as 
INTEGRAL. 

The radiation efficiency of an ESD (the fraction of the 
total gravitational energy lost to radiation) is low, tjesd = 
L^/r/oMc 2 w 1.6 x 10" 2 (m/50) _1 . The low efficiency im- 
plies that most of the photons are trapped and much of the 
gravitational energy dissipated by viscosity is advected into the 



BH. We call such a flow with super-Eddington accretion rate 
"photon-trapping accretion flow" (PTAF). 

The low conversion efficiency causes the hump luminosity 
to be weakly dependent on the accretion rate. Therefore, even 
with the uncertainties on /, the mass of the black hole can be 
estimated from the hump luminosity. Using eq. 14 and integrat- 
ing from the inner radius to the trapping radius we get 

M BH = 2.8 x 10 6 ( " L " . ) M , (30) 
\ 10 44 erg/s ) 

where vL v is the luminosity anywhere in the flat part of the 
hump, and we have used the approximation f* trap x 3 dx/ (e x — 
1) w J °° x 3 dx/(e x - 1) = tt 4 /15, x in = fw/kT(r in ), and 
£trap = hv/kT(r txap ). This method for determining the BH 
mass is unique since the the exact values of accretion rate and 
viscosity do not enter the mass determination. This method 
does not work for sub-critical accretion. 

It is worth pointing out that in the context of standard 
accretion disks, under the assumption that the viscous stress 
is proportional to radiation pressure, the radiation pressure- 
dominated region is thermally and viscous unstable because 
of the sensitive temperature dependence on viscosity. The 
optically thick ADAF has no such instability (under the a- 
description) because the surface density S is proportional to 
m (eqs. 11-13). 

4. Applications to NLS1s 

4. 1. The Eddington ratio in NLS1s 

There are several claims in the literature regarding the large 
L/L E dd in NLSls (e.g. Pounds, Done & Osborne 1995; Boiler, 
Brandt & Fink 1996; Laor et al 1997; Collin et al. 2002). 
However, the available mass estimates are remarkably few and, 
in most cases, highly uncertain. In particular, detailed reverber- 
ation mapping for the purpose of measuring the BH mass, are 
available for only a handful of sources. Several big campaigns, 
such as the one described in Collier et al. (2001) and Shemmer 
et al. (2001), provided uncertain mass estimates because the ob- 
served continuum variability was very small. Currently there 
are seven objects with FWHM(H/3)< 1500 km/sec and mea- 
sured BH mass. Their observed properties are listed in Kaspi 
et al. (2000) and discussed further in Peterson et al. (2000). 
According to Peterson et al. (2000), the BH mass in those 
sources is considerably smaller than the masses measured in 
comparable luminosity BLSls. This interpretation is still ques- 
tionable because of the uncertain BLR size used in those mass 
estimates. Collin & Hure (2001) and Collin et al. (2002) ex- 
amined accretion rates in AGN using the Kaspi et al. (2000) 
relationship. They find that half the objects in their sample are 
accreting close to the Eddington rate or at a super-Eddington 
rate. The NLSls in Collin's sample have the largest Eddington 
ratio. 

To further test this idea, we have used the compilation of 
Veron-Cetty et al. (2001) that contains a heterogeneous sam- 
ple of 54 NLSls with FWHM(H/3) < 2000 and B-band ab- 
solute magnitudes. The reverberation relation of Kaspi et al 
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L °l( £ Bol/ L Edd) 

Fig. 4. The Eddington ratio distribution in the Veron-Cetty et al. 
(2001) sample calculated as explained in the text. 



(2000) allows a determination of -Rblr as a function of lu- 
minosity over a large range in optical continuum luminosity. 
Following Wang & Lu (2001), we estimated BH masses in this 
sample using the 5 100 A monochromatic luminosity which we 
obtained by assuming a "typical" optical continuum of the form 
L v oc i/~ - 5 and (as in Kaspi et al. 2000) a bolometric lumi- 
nosity of L Bo \ ~ 9AL A (5100A). Given the measured FWHM, 
which is taken to be the one-dimensional gas velocity, we cal- 
culated BH masses assuming virialized BLR. Combining with 
the observed flux and redshift, we estimated the Eddington ra- 
tio for all 54 sources. The results are shown in Figure 4. 
Inspection of Fig. 4 shows that the mean of the sample is 



log 



V-^Edd/ 



= 0.08 ±0.31. 



(31) 



i.e, most of the NLS 1 are super-Eddington accretors (the two 
extreme cases are IRAS 05262+44 and NGC 4051). Moreover, 
the theoretical results presented earlier, as well as the results 
of Wang & Zhou (1999) and Ohsuga et al. (2002), suggest that 
in such sources the radiated luminosity is significantly lower 
than the accreted gravitational energy because of photon trap- 
ping. In particular, in PTAFs, to ^> L-q \/ -^Edd- It thus seems 
reasonable to assume that most NLS 1 s are super-Eddington ac- 
cretors. 

As shown in the previous section, super-Eddington accre- 
tion can be identified by the presence of a soft X-ray hump in 
the spectrum. The following section addresses the possible ob- 
servation of such features. 

4.2. Soft X-ray humps in NLS1s 

Soft X-ray humps in NLS 1 spectra have been discussed in ear- 
lier papers (Turner, George & Nandra, 1998; Leighly 1999, 
Vaughan et al. 1999, Comastri et al. 2001; Turner et al. 2001; 
Puchnarewicz et al. 2001; Collinge et al. 2001; O'Brien et al. 
2001; Boiler et al. 2001). There were several suggestions relat- 



ing such features to blends of X-ray emission lines. We con- 
sider those unlikely due to theoretical reasons (Turner, George 
& Netzer 1999) as well as the lack of direct evidence in re- 
cent high resolution XMM observations of some of the sources. 
Below we discuss the spectrum of RE J 1034+396 as an ex- 
treme case and address the possible connection to the theoreti- 
cal spectra shown in the previous sections. 

The ultra-soft X-ray radiation of RE J1034+396 (z = 
0.003) was first discovered by Puchnarewicz et al. (1995). The 
spectrum shown in Fig. 5 shows a clear soft X-ray hump ex- 
tending from log v = 16.5 (the lowest observed frequency) to 
17.2. The data were taken from Puchnarewicz et al. (2001). 

The fittings made by Puchnarewicz et al. (2001) to the spec- 
trum of RE J 1034+396 were based on the standard accretion 
disk model (Czerny & Elvis 1987). The resulting accretion rate 
is 0.3-0.7 M cr and the inclination 6 ohs = 60 - 75°. Such a 
high accretion rate is inconsistent with a thin disk model. For a 
standard accretion disk model, the maximum frequency of the 
spectrum is roughly proportional to to 1 / 4 and the luminosity 
is proportional to to. A strong soft X-ray hump in such disks 
requires very large to. This would give very large luminosity 
which is the reason why Puchnarewicz et al. (2001) had to as- 
sume high inclination (to decreases the total observed luminos- 
ity by a factor of cos 6* b s ). Even with such extreme values of 
to, there were difficulties fitting the observed spectrum. 

The ESD model provides a natural way to increase to with- 
out a proportional increase in the bolometric luminosity. To il- 
lustrate this, we assume a face on system with a = a c = 0.1. 
The mass of the central black hole can be obtained from equa- 
tion (29). Fig. 5 shows three such models fitted to the spec- 
trum of RE J1034+396 representing M BH = 2.25 x 1O 6 M , 
to = 5, 10, 20 and / = 1.1 x 10~ 2 , 5.5 x 10~ 3 , 2.75 x 10~ 3 . 
The maximum frequency of the hump is determined by r; n 
which, in this model, is fixed at r; n = 3. The minimum fre- 
quency is determined by the accretion rate, as seen from the 
theoretical spectrum in Fig. 5. However, in this case the low 
frequency part of the spectrum is not observed and thus the ac- 
cretion rate is unknown. There are only two parameters, Mbh 
and r; n , that determine the total flux in the hump. The hard X- 
ray spectrum (at frequencies larger than ^ max ) is determined by 
fM and a c . Since M is unknown, / is the real free parameter 
of the fit. 

We find that to = 5 cannot provide a satisfactory fit to the 
data. The combination of to = 10 and frh = 0.055 provides a 
good fit for the spectrum of and m = 20 is also consistent with 
the data given the lack of observations at low X-ray energies. 
The large observational uncertainties on the hard X-ray spec- 
trum translate to large uncertainty on the value of /. The best 
fit is M BH = 2.25 x 1O 6 M , fM = 0.055M cr , M = 10M cr 
and n n — 3.0. We have also computed a model with r in = 2.0 
and found very little change thus the dependence on r m over 
the fitted energy range is very weak. 

The BH mass obtained from the model can be compared 
with the mass obtained by using the Kaspi et al. (2000) relation- 
ship for i? B LR vs L 5 ioo, and the observed line widths. Using 
the data in Puchnarewicz et al. (1995) with FWHM(H/3)=1500 
Km/s, we get M = 4.1 x 10 6 M Q , in good agreement with the 
spectra-fit results. 
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Fig. 5. The soft hump and hard X-ray spectrum of RE J1034+396. 
The data are taken from Puchnarewicz et al. (2001). The dashed, 
solid and dotted lines represent ESD models with rh = 5, 10, and 
20 respectively, and frk = 0.03. In all cases r ln = 3 and Mbh = 
2.25 x 10 6 M©. The fits include Comptonization in both the ESD and 
the hot corona. 

The present model allows two different kinds of X-ray vari- 
ations. The hard X-ray spectrum is the Comptonized compo- 
nent of the soft hump flux. This means that intrinsic variations 
in the hump luminosity must result in correlated hard X-ray 
variations but not necessarily in optical-UV variations. On the 
other hand, the hard X-ray luminosity is very sensitive to the 
value of /. Time dependent variations in / will lead to varia- 
tions in the hard X-ray flux yet the soft hump canremain con- 
stant because of photon trapping effects. This can be tested ob- 
servationally. 

5. Discussion and Conclusions 

In this paper we explored the structure and the spectrum of ex- 
treme slim disks with hot coronae. We focused on the question 
of whether such systems can fit, simultaneously, the soft X-ray 
hump and the hard X-ray spectrum of NLSls. We presented 
a self-similar solution of ESDs with hot coronae and showed 
that their spectrum contains low energy humps produced by the 
ESD and a higher energy power-law continuum with a cutoff 
produced by the hot corona. We find that the structure of the 
hot corona, which we assumed to be Thomson-thin, is inde- 
pendent of the mass of the black hole and is fully specified by 
the parameters fm and a c . The ESD structure and spectrum 
(the "hump") is insensitive to the accretion rate, but sensitive 
to the mass of the black hole. This can explain the similarity of 
the hard X-ray spectrum of Galactic black hole candidates and 
NLSls despite of the very large difference in BH mass. Our 
present models apply torii> 2.5 and / < 0.0375 (m/10) _1 . 
The calculations do not include GR effect and apply only to 
non-rotating BHs. 

Slim disk fluctuations have been suggested by Mineshige 
et al. (2000) as the origin of the soft X-ray variations in such 
systems. The reason is the unusual energy balace 

^inag ^ Wgrav ^> Urad (32) 



where u denotes the energy density. The trapped photon ac- 
cretion flows discussed in this work are likely to show two 
different types of variations: 1) simultaneous variations of the 
soft hump and the hard X-ray without any significant time lag; 
2) small variability of the soft hump and larger variation in 
the hard X-ray spectrum. The first type is related to the so 
called "photon bubble instability" (Gammie 1999) that can be 
very important in SDs since the energy densities of the trapped 
photons and the magnetic field are larger than those in stan- 
dard accretion disks. In this case, the hard X-rays will follow 
closely the soft hump variations. Using the simplest version of 
Comptonization (Sunyaev & Titarchuk 1980) we find that the 
ratio of the hard X-ray flux to the soft hump flux is roughly 
p(p + 3)(2p + 3y 1 (vHx/vhum p y P where mx and f hump are 
representative frequencies of the hard and the soft (hump) X- 
ray flux. This ratio is sensitive to the hard X-ray spectrap index 
p. This type of variations agrees with at least some observations 
(e.g. Leighly 1999, Dewangan et al. 2002). Regarding the sec- 
ond type, here the fluctuations are due to changes in the accre- 
tion rate. Such fluctuations cause large hard X-ray variations 
since that radiation is sensitive to the accretion. On the other 
hand, the soft hump flux remains almost constant due to pho- 
ton trapping. This type of variations has been reported in the 
Chandra observations of 1H 0707-495 (Leighly et al. 2002). It 
will be interesting to apply a more detailed combined spectral 
and stability analysis of this type to the observed X-ray spec- 
trum of NLSls. 

The calculations presented here include the new ingredi- 
ents of a hot corona, yet they are still oversimplified in various 
ways. First, we assumed a continuous corona which is charac- 
terized by a single parameter /. However, patchy corona may 
be more reaslistic for SDs. Second, the hot corona gas is as- 
sumed to rotate with Keplerian velocity. This leads to an un- 
derestimate of the advection in the hot corona since the hot gas 
may be sub-Keplerian. This uncertainty can be aborbed into 
the factor /. A possible improvement is to use the optically 
thin advection-dominated self-similar solution (Narayan & Yi 
1994) as an approximation of the hot corona. We also note that 
for the case of standard accretion disk, the factor / approaches 
unity (Svensson & Zdzarski 1994). However, in that case feed- 
back to the cold disk can be neglectable since / < 1 in the 
super-Eddington case. Third, we do not take into account the 
irradiation of the disk due to the hot corona. Future work will 
have to treat the more general case of irradiation and try to 
solve, self-consistently, for the disk-corona interface e.g. via 
ionization instability (Zycki, Collin-Souffrin & Czerny 2000). 

Another limitation of the model is the neglect of GR effects. 
Thus, the treatment of the physical processes near the disk in- 
ner boundary, and the calculation of the emitted spectrum from 
this region, are not very accurate. As explained, global energy 
considerations show that the models presented here do not de- 
viate much from the more realistic cases since a large fraction 
of the radiation is emitted way outside of R in = 3, our stan- 
dard inner radius. In general we expect models with full GR 
treatment and similar m to show a very similar hump and to 
require a somewhat smaller / for the same hard X-ray flux. We 
will test these qualitative estimates in a future paper. Finally, 
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we neglected any spectral dependence on the inclination angle 
of the disk. This will be included in future models. 

Our new model enables a comparison with several already 
observed NLS Is. Study of the Veron-Cetty et al. (2001) sample 
showed that the mean Eddington ratio in this sample is larger 
than unity which suggests super-Eddington accretion rates in 
many NLS 1 . As an application of the model, we fitted the spec- 
trum of RE J 1034+396 and obtained a good agreement, within 
the observational uncertainties, over the entire observed X-ray 
band. We also find that the mass of the central black hole ob- 
tained from the model is in good agreement with the mass ob- 
tained by using the reverberation mapping method. 

Regarding new observations to test the model, perhaps the 
most important are Chandra and XMM spectra that cover the 
lowest energy range. This will enable a careful study of the 
"red" side of the ESD hump in several sources and a more di- 
rect determination of the accretion rate. A careful study of the 
"knee" spectral region is another observational challenge with 
implications to the type of black hole (via the value of J2j n ) and 
the energy production in the corona (via the value of /). 
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